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Abstract 



We calculate in the framework of two Higgs doublet models the differential 
decay rate for inclusive B — > r v X transitions to order 1 / m% in the heavy 
quark expansion, for both polarized and unpolarized tau leptons. In contrast 
to the situation in the standard model, we find sizeable corrections. A 

systematic heavy quark expansion significantly reduces the theoretical uncer- 
tainties in the calculation compared to the existing free quark decay model 
results. From the experimental measurement of the branching ratio, we derive 
the bound tan (3 < 0.51 mn± [GeV]. We point out that the tau polarization is 
potentially a more sensitive probe of multi scalar models than the branching 
ratio itself. 
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I. INTRODUCTION 



Semileptonic B meson decays into a tau lepton are sensitive to extensions of the standard 
model (SM) with several Higgs fields. In two Higgs doublet models (2HDM) [|I| the exper- 
imental measurement of the decay rate B — > r v X provides the strongest upper bound on 
tan/3/m#± for large values of tan/3. Theoretically, charged scalar contributions could 

increase BR (B — ► r v X) by as much as an order of magnitude compared to the SM 
This could have provided an exciting solution to the semileptonic B meson branching ratio 
problem []7|,§. While this scenario is ruled out by the recent ALEPH measurement 0, the 
ratio of decay rates into the tau and light lepton channels still provides a sensitive probe to 
multi scalar models. Moreover, the possibility to study the tau polarization offers a greater 
variety in both the experimental and the theoretical analyses Q. 

Recently, it has been observed that inclusive semileptonic decays of hadrons containing 
a single heavy quark allow for a systematic, QCD-based expansion in powers of 1/ttiq ||10|| . 
The heavy quark limit (mg — > oo) coincides with the free quark decay model and there are 



no corrections to this result at order 1 / tuq JT0|,[TT]] . The leading nonperturbative corrections 
are of order 1 / roq and depend on only two hadronic parameters, which parameterize forward 
matrix elements of local dimension-five operators in the heavy quark effective theory (HQET) 
1 372 1 - These corrections have been computed for a number of processes [ |IT| , |iT5| -|TE| . 

In the framework of the SM it was shown |18| that the heavy quark expansion provides 
more precise predictions than the spectator quark model for inclusive B — > r v X decays. 
The effect of including the nonperturbative 1/mg corrections is numerically small. More 
important is that HQET gives an unambiguous definition of the quark masses [l9j that 
determine the inclusive decay rates. Using the relation between the masses of the heavy 
quarks in terms of the HQET parameters leads to a significant reduction in the theoretical 
uncertainties [pj[| . 

In this paper, we use a combination of the operator product expansion and HQET to 
study the decay B — > tvX in the framework of 2HDM. In section 2 we define the model. 
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In section 3, we present the analytic results to order 1/m 2 for the branching ratio and for 
the tau polarization. In section 4, we study the numerical predictions for the total decay 
rate and for the tau polarization. Section 5 contains our conclusions. 



II. THE MODEL 

Extensions of the standard model scalar sector are constrained by the experimental value 
of p = m 2 ^ / (m| cos 2 9w) — 1 1 20 1 and by the strong limits on flavor changing neutral currents 



(FCNC) . The first of these constraints is naturally fulfilled if the Higgs sector contains only 
doublets (and singlets), providing the tree level prediction p — 1. FCNC are eliminated at 
tree level if we assume natural flavor conservation (NFC) ||21|| , a construction that requires 
fermions of a given electric charge to couple to a single Higgs doublet. The phenomenology 
of such models has been widely studied |L|. In these models weak decays are mediated 
besides the W bosons by charged scalar particles as well. 

For simplicity we consider 2HDM with NFC and comment on the more general case 
below. After spontaneous symmetry breaking only two physical charged scalars remain, 
H ± . Such a model has still several variants. If the Higgs doublet that couples to up-type 
quarks also couples to down-type quarks or to charged leptons then the experimental results 



on the decays B X s j |22| and Z — ► bb |23[ yield very strong constraints and the scalar 
contribution to B — > r v X is negligible. Therefore, we only consider models in which one 
Higgs doublet couples to down-type quarks and charged leptons and the other to up-type 
quarks (model II in the language of [l|]). The minimal supersymmetric standard model 
belongs to this class, so our results also apply to this model (up to radiative corrections). 

The Yukawa interaction of the physical charged scalars with fermions is determined 
by tan/? (the ratio of the vacuum expectation values of the two Higgs doublets), by the 
fermion masses and by the CKM matrix. The terms in the effective lagrangian relevant for 
B — > t v X decays are 

AG F 



C = -Vr 



do' 



V2 



(crfP L b) (f^P L is T ) - R (cP R b) (fP L u T ) , (2.1) 
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where 

R = r m T m b , r = , (2-2) 

m H ± 

and P^x, = ~(1 ± 75). We denoted by an upper index Y the running quark masses (this 
distinction is unimportant for the tau, see discussion in section 4). The first term gives the 
SM contribution, while the second one gives that of the charged scalars. We neglect a term 
proportional to m\: first, it is suppressed by the mass ratio mf/mj ; second, it cannot be 
enhanced by the possibly large factor tan 2 f3. 

Our analysis holds with minor modifications for general multi Higgs doublet models with 
NFC (for a recent comprehensive analysis, see In such models, instead of the single 
parameter tan/3, three complex coupling constants determine the Yukawa interactions of 
the charged scalars. The parameters X, Y and Z describe the couplings to up-type quarks, 
down-type quarks and charged leptons, respectively. Then in the effective lagrangian (|2.1|) 
r 2 = XZ*/m 2 H , where m# is the lightest charged scalar mass (assuming that the heavier 
charged scalars effectively decouple from the fermions). Thus, although r 2 is real in 2HDM, 
we treat it as complex so that our results apply to general multi Higgs doublet models. 



III. HEAVY QUARK EXPANSION 

In this section we present the analytic results necessary for our numerical analysis. The 
techniques by which they were obtained are described in detail elsewhere Hl~3|-|l^], so we give 
only a brief summary, followed by the results of the computation. 

The inclusive differential decay distribution is determined by the imaginary part of the 
time-ordered product of two flavor- changing currents, 

T = -ijd 4 xe- iq - x (B\T{O\x),O(0)}\B), O = c('fP L - mXrP R )b. (3.1) 

Since over most of the Dalitz plot the energy release is large (of order mi,), the time-ordered 
product can be written as an operator product expansion, in which higher-dimension oper- 
ators are suppressed by powers of A/m&, where A is a typical low energy scale of the strong 
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interactions. To this end it is necessary to separate the large part of the b quark momentum 
by writing p b = m b v + k, where v is the velocity of the decaying B meson, and its resid- 
ual momentum k is of order A. This separation is most conveniently performed by using 
the formalism of the heavy quark effective theory fll2| . One can then evaluate the matrix 
elements of the resulting tower of nonrenormalizable operators with the help of the heavy 
quark symmetries. 

The leading term in the expansion reproduces the result of the free quark decay model 
pT0|] , while giving an unambiguous definition of the heavy quark mass fL9|. The leading 



nonperturbative corrections are of relative order 1 / ml and can be written in terms of two 
parameters, Ai and A2, which are related to the kinetic energy K b of the b quark inside the 
B meson, and to the mass splitting between B and B* mesons: 

K b = -^-, m%. -m| = 4A 2 . (3.2) 

The operator product expansion for the SM contribution has been presented in refs. 
I7|,IB]. Including the charged Higgs contribution is, in principle, a straightforward general- 



ization. We shall present only the final results. The tau lepton can have spin up (s = +) 
or spin down (s = — ) relative to the direction of its momentum, and it is convenient to 
decompose the corresponding decay rates as 

r (B -> t(s = ±) uX) =-r±f. (3.3) 

The total rate, summed over the tau polarizations, is given by T, while the tau polarization 
is A po \ = 2r/r. Let us further decompose 

\V cb \ 2 G%ml 



r w + l -^-r H -2Re(R)^r I 
4 



(3.4) 



1927T 3 

and similarly for T. The subindices W, H, and / denote the W mediated (standard model), 
Higgs mediated and interference contributions, respectively. This notation is particularly 
convenient, because in the spectator model IV = r#, Th = —Tw and Ti = 0. In the 
m T , m c — > limit of the spectator model Tw = = Tj = 2Tfj = — 2T\y = 1. After 



integration over the invariant mass of the lepton pair and the neutrino energy, the differential 
decay rate depends only on a single kinematic variable E T , which denotes the tau energy in 
the rest frame of the decaying B meson. We will use the dimensionless variables 



2E^ 



y 



rrr, 

p= — b 

mi 



mi 



Pr 



mi, ml m b 

The standard model contribution to the lepton spectrum (dT w /dy and dT w /dy) has 
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(3.5) 



been calculated in ref. ||18|| , and we do not present it here. In terms of these quantities, 
however, the Higgs mediated contribution is quite simple. The reason is that the kinetic 
operator does not violate the heavy quark spin symmetry, thus the corrections proportional 
to Ai are the same for the W and Higgs mediated terms. We obtain 

dT H 



dy 

df H 
dy 



dT w A 2 24x (l - x )[y(3-y) - p r (8 - 3y)] x /y 2 - Ap T 
dy m 

dT\y 



{ l + p T -y 

A 2 12x (l - x Q ) (y 2 - 4p r )(3 -y-p T ) 



dy 



mi 



where 



x = 1 - 



1 + Pt - y 



For the interference term we get: 



dVi 
dy 



v /y 2 -Ap T {6x 2 (l + p T -y) 
Ai 



m 2 b {l + p T -y) 



3A 2 

m 2 b (l + pT -y) 



and 



dy 



2(y 2 - 4 Pr ) + 2x {y 2 - 3j/(l + p T ) + 8p T ) 
-4 (3 + 14p T + 3p 2 T - 9y(l + p T ) + Ay 2 ) 
2x (12 + 8p r - 3y{5 + p T ) + Ay 2 ) 
-xl (15 + 10p r + 3p 2 T - 3y(7 + 3p T ) + 8y 2 ) ] } 

A 2 3x (l-x ) (y 1 - 4p T )(2 - y) 



mi 



l + pr - y 



(3.6) 
(3.7) 



(3.8) 



(3.9) 



(3.10) 



The explicit form of the total decay rates are obtained by integrating the above expres- 
sions over the kinematic range 2^fp~^ < y < 1 + p T — p. We do not reproduce here the explicit 
expressions for the standard model contributions T\y and Fyy that can be found in ref. |jTB| . 
For the Higgs mediated term we obtain: 
12A 



T H = T W + ^{VX [2(1 - p T ) 2 + p(5 - p + 5p r )] - 6p(l - p 2 T )A - 6p(l + p 2 )M , (3.11) 

777.L I- J 
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f H = -t w + —j\ [(1 - m T ) 2 - p] 2(1 - p r )(l + m T f - p 



m 



+p(5 — 3m T — p r + 15p r m T ) /(l — m T ) 
-2p(3 + 2p T -4p r p-5p 2 )ln- 



2,!„ (1 ~^r) 2 



(3.12) 



where we defined 
A = In 



2m T p 



B = In 



2m T 



(3.13) 



l-p T ) 2 -p(l + p T )-(l-p T )v / A' " l + p T -p+VA 
A = 1 — 2(p T + p) + (p r — p) 2 , and rh T = m T /mb = y/p\. For the contribution of the 
interference term we obtain 
Ai 



r,= Va{(i 



2m 2 



l-5p-2p^ + 10p r -5p T p + p^ 



3A 9 



2m 2 



11 - p + 2p T + 8p 2 + 35p r p - p 2 )} 



+6 1 + 



Ai 
2m 2 



p 2 (l - p T )A + (p 2 + 2p T )(l + p T ) - Ap T p)B 



9A 5 



p(2 + p - 5p T p - 2p 2 T )A + 2p(l - p T f - 6p T + 2p 2 + p 2 + 5p T p 2 £ , (3.14) 



and 



3A 9 



2m 2 



[(1 - m T f -p] (1 - p T )(l + m T f - p 2 - 2m T p(l - 5p T )/(l - m T ) 



-2p(l-p-p T )(l + 3p T )ln 



1 ~ rhr) 
P 



(3.15) 



In the limit Ai,A2 — > 0, corresponding to the free quark decay model, our results agree 
with |J. However, we disagree with jj| on the size of the interference term. There it was 
presumably multiplied with an extra factor of m T jm} ) . As expected, the parameter Ai only 



enters the total decay rate proportional to the spectator model result, namely the total 
decay rate depends on Ai only through the combination (1 + Xi/2m b i ). By virtue of the 



coupling of the charged scalar to the lepton pair (|2.1|), the Higgs mediated term produces 
tau leptons with positive polarization. This indicates that, in addition to the decay rate, 
the tau polarization is also sensitive to new physics. 

IV. NUMERICAL RESULTS 

1. Input parameters 

When we neglect the tiny contribution from b — > u transitions, the parameters entering 



our calculations are the mass of the tau lepton m T = 1.777 GeV [24], the heavy quark masses 



mi, and m c , the hadronic parameters Ai and A2, the quark mixing parameter \V c b\, and the 



combination of the unknown Higgs sector parameters R (2.2). To utilize the full power of 
the heavy quark expansion it is important that not all of these parameters are independent 
18|,25]. In fact, the same HQET parameters Ai and A2 appear in the expansion of the heavy 



meson masses in terms of the charm and bottom quark masses (up to the running of A2) 



1 t A! + 3A 2 

m B = m b + A h . . . , 

2m b 

t Ai -|- 3Ao , x 

m D = m c + A - 1 + . . . . (4.1) 

The parameter A can be associated with the effective mass of the light degrees of freedom 
inside the heavy meson [^,|T^][|. For each set of values {A, Ai,A 2 }, eq. ( |4.1| ) determine m c 



and nib- The consistency of the heavy quark expansion requires that these values of the 



*It has been argued recently [28] that the definition of A is ambiguous. If this ambiguity turns 
out to be numerically significant that could probably be accounted for by taking a larger range for 
A in our numerical calculations and would slightly increase the uncertainties. 
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quark masses are used in the theoretical expressions for the decay rates. Using only three 
independent parameters {A, Ai, A 2 } instead of the four quantities {m&, m c , Ai, A2} reduces 
the theoretical uncertainties significantly. 

These HQET parameters are genuinely nonperturbative. While the value of A2 is related 
to the B* — B mass splitting A2 = (jn 2 B* — ^b)/4 — 0.12 GeV 2 , there is no similarly simple 
way to determine A and Ai. We expect this value of A2 to be accurate to within 10%, as a 
result of the finite b quark mass and the experimental uncertainties. The parameters A and 
Ai can only be estimated in models of QCD at present, or constrained from phenomenology 



25fl . QCD sum rules have been used to compute Ai |29|,|3(|, but these calculations suffer 
from large uncertainties. There is increasing theoretical evidence that, in accordance to its 
definition ( |3.2| ), Ai is negative [pi], and its magnitude cannot be too large |32|. Here we use 



< — Ai < 0.3 GeV , which is also supported by a recent QCD sum rule calculation |33 
Assuming Ai < 0, the phenomenological analysis of ref. |25| prefers values of A lower than 
QCD sum rules f2||g. Here we take 0.4 < A < 0.6 GeV. 

Finally, a subtle point is that of the perturbative QCD corrections. While we know 
exactly the 0(a s ) corrections to the W mediated term in the spectator model, analogous 
calculations for the H and I terms have not been carried out. For the W term t] T ~ 0.90 |B5 



we shall also use the 0(a s ) correction to B — > evX : i] e ~ 0.88. The leading corrections to 
the Higgs coupling are incorporated by running the Yukawa coupling. Therefore, R depends 
on the running mass of the b quark mft(/z). We take a conservative range 0.9 rrif, < mj < 
1.05 mi, corresponding to m c < fi < m^. For the interference term we include an extra 
(—5 ± 10) % correction to represent the uncertainty in the QCD correction. 

To conclude the above discussion we summarize the ranges for the various input param- 
eters that we use in our analysis: 

0.4 < A < 0.6 GeV, < -Ai < 0.3 GeV 2 , 

(4.2) 

0.9m fe <ml < 1.05m b , 0.11 < \ 2 {m b ) < 0.13 GeV 2 . 
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2. Branching ratio 



Normalizing the branching ratio in the tau channel to that into light leptons separates the 
theoretical and experimental uncertainties in the numerical analysis, eliminates the otherwise 
significant uncertainties from the values of \V c b\ 2 and m| , and also reduces the sensitivity to 
the unknown QCD corrections. We plot the theoretical prediction for the branching ratio 
as a function of r in Fig. 1. Our result is given by the shaded region between the solid lines. 

Using the recent measurement of the branching fraction into final states with a tau lepton 
HI and that with an electron 1 12011, 



BR{B ^tvX) = 2.76 ±0.63%, 

BR(B -> evX) = 10.7 ±0.5%, (4.3) 
we obtain the la upper bound 

r<0.51GeV" 1 . (4.4) 

The 95% CL upper bound is r < 0.55 GeV -1 . 

There are a number of points to be made regarding this result: 

a. The 1/ml corrections to the free quark decay model turn out to be more significant 
than in the SM: while the W mediated contribution is suppressed by —4% and —8% in the 
electron and tau channels respectively, the Higgs mediated term is enhanced by approxi- 
mately 14%, and the interference term is suppressed by as much as 25%. As a result, the 
difference between the 2HDM and the SM becomes more significant than expected from the 
free quark decay model. 

b. The uncertainty in our result (corresponding to the width of the shaded region in 
Fig. 1) takes into account the theoretical uncertainties in ( [4.2| ), as well as an estimate of 
the 1/rriQ and more importantly the perturbative QCD corrections. If the latter terms were 
known, the width of the shaded region would be reduced by about a factor of two for not 
very small values of r. 
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c. The improvement of our calculation over the spectator model originates from large 
1 / m 2 , corrections, and from using mj, and m c as determined by eq. ( |4.1| ) rather than treating 
them as independent input parameters. To illustrate this, we plotted with dashed lines in 
Fig. 1 the prediction of the free quark decay model corresponding to 1.4 < m c < 1.5 GeV 
and 4.6 < < 5 GeV. Our predictions are not very sensitive to the values of A and Ai. 
Adopting the range 0.3 < A < 0.7 GeV and —0.5 < Ai < 0.5 GeV 2 , which covers a range 
for m c and m& larger than the ones previously mentioned, makes the allowed band slightly 
wider than in Fig. 1. 

d. As the interference term gives rise to a dip in the branching ratio for small values of r, 
this measurement cannot probe values of r smaller than 0.43 GeV" 1 . Within the spectator 
model, this limit is r < 0.6 GeV -1 . 

e. In addition to our improvement of the theoretical calculation we would like to point 
out that some of the earlier analyses quoted too strong bounds on r by (a) neglecting the 
uncertainties in the quark masses; (b) underestimating the size of the interference term 
by about a factor of three based on and (c) overestimating the SM prediction for the 
branching ratio. 

3. tau polarization 

The polarization of the tau lepton, A po \ = 2T/T, being a ratio of decay rates, is subject 
to much smaller uncertainties than the rates themselves. We find that the numerical value 
for A P oi is rather insensitive to variations in A and Ai. Allowing these parameters to vary 
within the ranges ( |4.2|) , we find the allowed values of the tau polarization as a function of 
r. Our result is given by the shaded region between the solid lines in Fig. 2. 

A few points are in order regarding this result: 

a. The 1/mf corrections to the free quark decay model increase ^4 po i for small values of 
r. More importantly, these corrections make the tau polarization almost a monotonically 
increasing function of r. While the interference term contribution to T is small (it is exactly 
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zero in the spectator model), it is important that the contribution of the interference term 
to T is significantly suppressed while that of the Higgs mediated term is enhanced due to 
the 1/ml corrections. 

b. The uncertainty in our result (corresponding to the width of the shaded region in 
Fig. 2) takes into account the same theoretical uncertainties as in Fig. 1. The most important 
of these are the perturbative QCD corrections: if they were known, that would make the 
prediction for A po \ accurate at the 1% level for r < 0.3 GeV -1 . 

c. Besides the suppression of the interference term and the enhancement of the Higgs 
contribution to the total rate, the improvement of our calculation over the spectator model 
is again due to using rrif, and m c as determined by eq. ( |4.1| ) rather than treating them as 
independent input parameters. For comparison, we plot (with dashed lines) the prediction 
of the free quark decay model corresponding to the same ranges of the quark masses as in 
Fig. 1. 

d. Most important is that a tau polarization measurement of an accuracy of about 10% 
would be sensitive to values of r as small as 0.3 GeV -1 . As pointed out in the previous 
section, such a small value of r cannot be probed with the total branching ratio. The 
theoretical prediction for A po \ is accurate at the 1% level in the standard model —0.70 < 

< —0.71 jl8| . A precise measurement of this quantity would be a very stringent test 
of the SM, or alternatively provide a determination of the parameter r. 

Finally we would like to compare the bounds that can be achieved from the inclusive 



decay B — > r v X to that from the exclusive decays B — > £ v. In ref. |36[ it was shown that the 
purely leptonic decays of B mesons are enhanced compared to the SM for r > 0.27 GeV -1 , 
and a bound r < 0.52 GeV" 1 has also been claimed. It is important to note that deriving a 
bound on r from this decay requires an independent estimate of \V u i\ and the B meson decay 
constant Uncertainties related to the values of these parameters have been neglected 



in 



36| . These uncertainties are, however, very significant; taking them into account (fs > 
140 MeV and \V u b\ > 0.0024) allows one to obtain only r < 0.80 GeV -1 . Concerning future 
bounds from this process we estimate that an enhancement of BR (B — > £ v) will only be 
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observable if r > 0.35 GeV -1 . We conclude that a stronger bound on r will be provided by 
a measurement of the tau polarization in the inclusive B — > tu X decay, while at present 
the best bound comes from the measurement of the branching ratio in the same decay. 

V. SUMMARY 

We investigated in detail the effects of a charged scalar on the inclusive decay B — > tv X. 
We used the heavy quark expansion to incorporate the nonperturbative 1/rajj corrections, 
and the relation between the heavy quark masses implied by the HQET. This resulted in 
a reduction of the theoretical uncertainties by more than a factor of two compared to the 
calculation in the free quark decay model. While the I /ml corrections are usually less 
than 5-10% in standard model calculations, in the present case they turned out to be more 
significant. To derive a bound on r, special care has to be taken of the uncertainties in 
the theoretical calculation which have been mostly neglected so far. From the presently 
available experimental data the strongest bound that can be obtained is r = tan (5 / m H ± < 
0.51 GeV -1 . We pointed out that the tau polarization is more sensitive to small values 
of r than the branching ratio. While BR(B — > tvX) is insensitive to r < 0.4 GeV -1 , 
a measurement of A po \ could probe r ~ 0.3 GeV -1 or even lower. This is particularly 
interesting since the recent ALEPH measurement of the branching ratio excludes very large 
deviations from the standard model. We hope that this will encourage experimentalists to 
measure the tau polarization in B decays with the best achievable accuracy. 
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FIGURES 



FIG. 1. r (B — > tuX)/T{B — > ePl) as a function of r = tan.j3/m H ±. The shaded area 
between the solid lines is our result. The area between the dashed lines gives the free quark decay 
model result. The dash-dotted lines give the experimental la bounds. 

FIG. 2. The tau polarization, A po \ = 2T/T, as a function of r = tan (3/mjj± ■ The shaded area 
between the solid lines is our result. The area between the dashed lines gives the free quark decay 
model result. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403376vl 
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